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In this paper, we derive comparison results for terminal values of 
d- dimensional special semimartingales and also for finite-dimensional 
distributions of multivariate Levy processes. The comparison is with 
respect to nondecreasing, (increasing) convex, (increasing) direction- 
ally convex and (increasing) supermodular functions. We use three 
different approaches. In the first approach, we give sufficient con- 
ditions on the local predictable characteristics that imply ordering 
of terminal values of semimartingales. This generalizes some recent 
convex comparison results of exponential models in [Math. Finance 8 
(1998) 93-126, Finance Stock. 4 (2000) 209-222, Proc. Steklov Inst. 
Math. 237 (2002) 73-113, Finance Stoch. 10 (2006) 222-249]. In the 
second part, we give comparison results for finite-dimensional dis- 
tributions of Levy processes with infinite Levy measure. In the first 
step, we derive a comparison result for Markov processes based on a 
monotone separating transition kernel. By a coupling argument, we 
get an application to the comparison of compound Poisson processes. 
These comparisons are then extended by an approximation argument 
to the ordering of Levy processes with infinite Levy measure. The 
third approach is based on mixing representations which are known 
for several relevant distribution classes. We discuss this approach in 
detail for the comparison of generalized hyperbolic distributions and 
for normal inverse Gaussian processes. 

1. Introduction. Stochastic ordering and comparison results for stochas- 
tic models are topics of particular interest which have undergone an intensive 
development. The main aim is to derive approximation results and bounds 
for complex stochastic models where fundamental model assumptions are 
only partially known. Important fields of application include insurance math- 
ematics, risk management, reliability theory, renewal theory, stochastic net- 
works, financial mathematics, statistical physics and many others. Several 



Received October 2004; revised June 2005. 

AMS 2000 subject classifications. Primary 60E15; secondary 60G51, 60G44. 
Key words and phrases. Convex ordering, jump diffusion process, Levy processes, prop- 
agation of convexity. 

This is an electronic reprint of the original article published by the 

Institute of Mathematical Statistics in The Annals of Probability. 

2007, Vol. 35, No. 1, 228-254. This reprint differs from the original in pagination 

and typographic detail. 



1 



2 



J. BERGENTHUM AND L. RUSCHENDORF 



comprehensive books survey wide areas of this research (see [15, 17, 22, 23, 
24]). Most of the development is concentrated on one- or multi-dimensional 
orders. There is also an extended theory for ordering of discrete-time pro- 
cesses (like queuing sequences, renewal sequences, Markov chains) or related 
point processes with a wide variety of applications. Stochastic ordering re- 
sults for the classical stochastic order (w.r.t. nondecreasing functions) have 
also been established under various conditions for diffusion-type processes 
(cf. [6, 7, 11, 18, 25]) and for Markov processes (cf. [16]). These results are 
parallel to classical comparison theorems for solutions of differential equa- 
tions. Some interesting convex comparison results for exponential stochastic 
models have been developed in recent financial mathematics papers (see 
[2, 3, 5, 8, 9, 10]). The main aim in these papers is to derive sharp upper 
or lower bounds for option prices. The methods used in these papers are 
based on stochastic calculus (Ito's formula) and the propagation of convex- 
ity property (see [2, 3, 5, 8]) as well as on the coupling method (see [9, 10]). 

In our paper, we first give an extension of the approach in [2, 3, 5, 8] 
to obtain comparison results for multidimensional semimartingales for var- 
ious diffusion and monotonicity stochastic orders. This leads to stochastic 
ordering results of the following type. Suitable ordering of the local char- 
acteristics plus a "propagation of ordering" property imply that the basic 
semimartingale S can be compared to a Markovian semimartingale S* at 
any fixed time T. This result is worked out in Section 2. The main tool for 
the proof of the comparison result is the supermartingale property of the 
comparison process Q(t, St) based on the propagation operator Q(t, s), which 
for an ordering function g, is given by 

g(t,s) = E*(g(S* T )\S; = s). 

Our comparison result has several applications. We state several comparison 
results between diffusion processes, stochastic volatility models, diffusions 
with jumps and Levy processes. As a very specific application, we obtain a 
complete list of the ordering properties of multivariate normal distributions 
w.r.t. diffusion and monotonicity orders. 

In Section 3, we use a different technique to obtain ordering results for 
the finite-dimensional distributions of Levy processes. For the comparison, 
we use an approximation result to reduce the ordering problem to the com- 
parison of compound Poisson processes. These can be ordered by a direct 
coupling argument and by a general ordering result for time-homogeneous 
Markov processes based on the existence of a monotone separating transition 
kernel. The ordering conditions are then formulated in terms of the Levy 
measures as well as the drift and diffusion coefficients. As an example, we 
discuss ordering properties of normal inverse Gaussian (NIG) processes. 

Finally, in Section 4, we derive ordering results based on mixing-type rep- 
resentations which are known for several interesting classes of distributions 
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like generalized hyperbolic (GH), multivariate t and elliptically contoured 
distributions (cf. [4] for applications to risk management). We discuss this 
method in detail for the particular case of GH-distributions. Our main re- 
sults describe the ordering properties of GH-distributions in terms of the 
parameters of the mixing variables, which are identical with the basic pa- 
rameters of the distributions. For NIG-processes, these results also imply 
ordering results for the finite-dimensional distributions. 

2. Comparison results for semimartingales. We derive comparisons of 
terminal values of multivariate special semimartingales with absolutely con- 
tinuous characteristics and also of finite-dimensional distributions of Levy 
processes. The comparison is based on stochastic orderings <jf induced by 
the following function classes F: 



F st 

F :x 

(1) -^dcx 
F 3m 
•^icx 



{/ 
{/ 
{/ 
{/ 



i, f is increasing}, 
I, f is convex}, 
I, f is directionally convex}, 
I, f is supermodular}, 



•Ft F jt, -^idcx • — F i cx n F jt, F[ sm . — F 3m n F s 



or by subclasses of them, in particular by single elements g of these classes. 
We consider the componentwise order on M. d , that is, for x,y £ M. d , x < y 
holds if x % < y % for all i < d. For random vectors X and Y, the stochastic 
ordering X < T Y is defined by Ef(X) < Ef(Y) for all f € F such that 
f(X) and f(Y) are integrable. For the specific classes F in (1), we denote the 

orderings X <jr Y by X < st Y, X < cx Y, X <d cx Y, Orderings induced 

by these function classes are also generated by all functions from T n C°° , 
where F denotes one of the classes in (1) and C°° is the set of infinitely 
differentiable functions, as well as by many other order generating function 
classes Fq C F. For results on stochastic orders, we refer to [15, 17, 22] and 
[24]. 

For a finite time horizon T < oo and d E N, we consider (i-dimensional 
special semimartingales S and S* on stochastic bases (Q,A, (At)t£[o,T]i P) 
and (fl*,A*, (A.f)te[o,T])F'*), respectively, with predictable characteristics 
(B(h),C,v) and (B*(h*), C*, u*), respectively, where the drift components 
B(h), B* (/i*) depend on the choice of truncation functions h, h* . We assume 
throughout the paper that the characteristics are absolutely continuous and 
denote the corresponding differential characteristics by b(h),c and K, where 
b(h) is a predictable d-dimensional process, c = (c t3 )ij<d is predictable with 
values in the set of symmetric, positive semidefinite d x d matrices with 
real entries M + (d, M) and K^jidx) is a transition kernel from (0 x [0,T],"P) 
into (R d ,B d ) that satisfies K W)t (Q)) = and / K Ujt (dx)(\x\ 2 A 1) < 1 (cf. [13], 
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Proposition II. 2. 9). For A,B G M+(d,M), the positive semidefinite order 
A < psd B is defined by x T (B - A)x > for all x G R d . 
By 

S~(b(h),c,K)~(b(h),c,K) h 

we indicate that 5 1 has predictable differential characteristics b(h),c and 
if. Special semimartingales are characterized by (\x\ 2 A |x|) * z/ G A\ oc and 
/i = id may serve as "truncation function" (cf. [13], Proposition II. 2. 29). We 
define b := 6(id). For example, Levy processes with existing first moments are 
special semimartingales with deterministic differential characteristics that 
are given by 

(2) b t (uj) = tb, ct(w)=i£, K t (u)=tF, 

where b G R d and £ G M(d,M) are constant and F denotes the Levy measure 
(cf. [13], Corollary II. 4. 19, [21], Corollary 5.25.8). In this case, we write 
S ~ (6, Ti,F). Another example for a special semimartingale is a diffusion 
with jumps that is given by an SDE 

dS t = b{t, S t ) dt + a{t, S t ) dW t + <f>(t, S t - , y)(N(dt, dy) - q(dt, dy)), 

So = x, 

where W is a d-dimensional Brownian motion and iV is a Poisson random 
measure on M + x R d with intensity measure q(dt,dy) = dt (g) X(dy), that 
satisfies (|x| 2 A|x|)*z/G A\ oc . In this case, S has characteristics 5 1 ~ (6, c, if), 
with 

(3) c = aa T , K t (s,A) = jl Am (<P(t,s,y))X(dy)=:X^ s '-\A); 

see [13], Remark III. 2. 29. 

For the special semimartingale 5*, we additionally assume a Markovian 
structure in the sense that its differential characteristics are deterministic 
functions of time t and S*_, that is, 

S* ~ (b* (t , 5;_ ) , c* (t , St ),K*(t, st ) ) . 

For notation and results on semimartingales, we refer to [13]. 

The first approach relies on an integro-partial-differential equation for the 
propagation operator Q(t,s). If Q(t,S*) is a local (^4^)-martingale under P* 
for some Q G C l,2 ([0, T] x R rf ), then under an additional condition, it will 
be established in the following lemma that the propagation operator Q(t,s) 
satisfies the integro-partial-differential equation 

D t 0(t, a) +£)D i 0(t, s) 

i<d 

(4) 

+ \ E VlMt,s)c*v(t,s) + J d (Ag)(t,s,y)Kts,dy)=0, 
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where (A0)(t, s,y) := G(t, s + y) - Q(t, s) - EiKd^iO^ s)y\ yeR d and f 
denotes the ith component. For the jump measure \i of S and 

(5) W(u,t,y):=Ag(t,S t -(u;),y), 
the integral process (W * fi)t is defined by 

(W*fJ,) t (u) := / W(uj,u,y)fj,(u;;du,dy) 

J[0,t]xR d 

= / Ag(u,S u -(u;),y)[i(uj;du,dy). 

J[0,t]xR d 

The predictable function W* and the integral process (W* * jjL*)t are defined 
similarly. The proof of the following lemma is given in the Appendix: 

Lemma 2.1 (Kolmogorov backward equation). Let S± ~ (b*(t,Sf_), 
c*(t,S*_),K*(t,St-)) be a d- dimensional special semimartingale, let 
g € C 1,2 ([0,T] x M d ) and assume that G(t,S^) is a local (Af) -martingale 
under P* . If (\W*\ * fi*)t G A^~ oc or if Q(t, •) is convex, then G(t,s) satisfies 
the Kolmogorov backward equation (4). 

In the sequel, for g G J 7 , we consider the functional G(t,s) defined by the 
propagation operator 

(6) g(t,s) = E*(g(S* T )\S; = s), 

where E* denotes the expectation with respect to P* . A crucial assumption 
in the comparison results for special semimartingales is the propagation of 
order property: for some ordering function g = G(T,-) G JF of interest, we 
assume that the ordering at time T is propagated to earlier times t £ [0,T], 
that is, g(t, •) G T . This order-propagating property is known in exponential 
diffusion models in financial context for T = JF CX as "propagation of convex- 
ity." The no-arbitrage price process 0(t,s) of a European contingent claim 
with convex payoff function g(s) = g{T,s) is convex in s, for all t G [0, T] 
(cf. [2, 3, 5, 8, 10]). 

Assumption PO(g) (Propagation of order). For some g G we say 
that S* satisfies the propagation of order property PO(g) if g(t, •) G T for 
< t < T. Similarly, S* satisfies PO(J zr ) for some JF C T if PO(s) holds for 
all g G Tq. 

In the following, we need some integrability properties of G(t,St). We 
define the classes 7 = T(S, S*) and T = F(S, S*) by 

-f:={g£F:g(t,s) GC^QO.T] x R d ), 
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Q(t,St) bounded below and EQ(t,St) < oo Vi G [0,T]}, 
F := {g E F -.g(t,s) E C^i^T) xR d ) and 

Q(t,St) is a process of class (D)}. 

If, for example, (7 has bounded linear growth, \g(x)\ < K\x\ for K < 00 
and 5, 5* are Levy processes or exponential Levy processes, then the inte- 
grability condition EG(t, S t ) = J J g(s + s*)P* s r-t (ds*)P St (ds) < 00 follows 
from the integrability of St and S^_ t . 

The general comparison result is formulated in the following theorem for 
the increasing directionally convex and for the increasing convex orders gen- 
erated by ^ r = ^ r idcx and T = T lCK . For the other order- generating function 
classes in (1), similar results hold true (see Remark 2.3). The proof of the 
following theorem is based on ideas similar to those in [2, 3, 5] and [8], which 
are concerned with stochastic exponential models. 

Theorem 2.2 [Increasing (directionally) convex comparison of semimartin- 
gales] . Let S, S* G Sp and Sq = Sq have differential local characteristics 
5~ (b,c,K) id and S* ~ (b*,c*,K*) id . Let W be as in (5) and let W* be 
defined similarly. 

1. Let g G Fidcx or g € J-idaa let Assumption PO(g) hold and assume that 
the associated processes (\W\ * ff)t and (\W*\ * fi*)t are in A^. Then the 
comparison of the differential characteristics 

(7) 6j( w ) < b*\t,S t .(co)), c\ j {u) < c"^,S t _(u,)), 

(8) / f(t,St-{u),x)K u , t {dx)< I f(t,S t .(u J ),x)K;(S^(u),dx), 
Jm d JR d 

\ x P-a.e. for all f G R + xl^xl^-tl with f(t, s, •) G J^dcx such that the 
integrals exist, implies that 

Eg(S T )<E*g(S* T ). 

2. Let g G J-i cx and let Assumption PO(g) hold true. Then the comparison 
of the differential characteristics 

(9) bi(u) <b*% s t _H), (^M^^psdCc^'^St-H))^, 

(10) / f(t,St-(u>),x)K Uit (dx)< [ f(t,St-(u),x)Kj(St-(u),dx), 

JR d JW. d 

\ x P-a.e., for all/el+x^xl^l with f{t,s,-) eT cx , such that the 
integrals exist, implies that 

Eg{S T ) < E*g(S* T ). 
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Proof. We prove that the comparison process Q(t, St) is an (y4t)-super- 
martingale under P. Then using the fact that Sq = Sq, we obtain the stated 
comparison result 

Eg(S T ) = EQ(T, St) < 0(0, So) = E*g(S* T ). 

1. Let g G Tide* or J^idcx- As (|W| * fi)t G A^ c , Ito's formula implies, sim- 
ilarly as in the proof of Lemma 2.1, that 

G(t,S t )=G(0,S )+M t 

+ [ {B t g( u ,s u .) + Y,^G(u,s u .)bi 

+ \ V Df,g(«,5 u _)c^+ / (A^s^,^^)}^, 
^ J 

where Mt := X^<dJ[o t] >?«-) gL/V^ is a local („4j)-martingale under P 

and JV£ denotes the ith component of the martingale part of the canonical 
decomposition of St = So + N t + B t . As Q(t, s) satisfies the integro-partial- 
differential equation (4), we obtain 

g(t,s t ) =g(o,s ) + M t 



(11) 



+ 2 E D^(u,S n _)(c^-c">,S u _)) 



+ / (Ag)(n,S M _,rr)(A: u (dx)-A::(5 u _,(ix)) 

JR d J 

=: G(0, S ) + M t + V t , 



where Vt is a predictable process of finite variation and, therefore, Vt G -4i oc - 
We now prove that inequalities (7) and (8) imply that Q(t,St) is a local 
(^li)-supermartingale under P. By Assumption PO(<?), Q(t, ■) is increasing 
and directionally convex for all t G [0,T]. Therefore, T>iG(t,s) and D?-f?(t,s) 
are nonnegative and inequalities (7) imply that the first and second terms 
of the integrand of Vt are nonpositive. For fixed (uj,u) G (fi, [0,T]), we define 
the function 



(12) 



T(x) :=(A^)(«,^_H,x). 
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As DfjT(x) = DfjQ(u, S u _ + x) > 0, for all x, T is directionally convex and, 
therefore, (8) implies that the last term of the integrand of Vt is nonpos- 
itive. This yields — Vt £ A^~ oc and it follows that Q(t,S t ) is a local (At)- 
supermartingale under P. 

It remains to prove that Q{t, St) is an (>l i )-supermartingale. If g € J^idcx is 
bounded below, Mt is bounded below and, therefore, is an (_4.t)-supermartingale 
under P. It follows that Q{t,St) is a super martingale as it is integrable. In 
the case where g £ Jidcx, we consider a localizing sequence r n for Q(t,St). 
Since for all t£ [0, T], we have P-a.s. that (G(t,St)) Tn — > Q{t, St), n— > oo, 
and Q(t, St) is of class (D), the convergence takes place in L 1 and, therefore, 
Q(t,St) is an („4 t )-supermartingale under P. 

2. Let g 6 ^"icx- Similarly to the proof of the first part, we show that 
Q(t, St) is an (^)-supermartingale under P. The evolution of Q{t, St) under 
P is given in (11). Observe that as in Lemma 2.1, convexity of Q{t, ■) implies 
(\W\ * /j)t,(\W*\ * fi*)t G A-ioc- As G(t,-) is increasing, comparison of the 
drift characteristics in (9) implies that the first term of the integrand of 
Vt is nonpositive. Due to the positive semidefiniteness of (c* lJ — c lJ )ij<d := 
(c* tJ (u, S u -)(uj) — c^(io))ij<d for fixed (u>,u), its spectral decomposition is 
given by (J2k<d -^fe e fc e i)*J<rfi where the eigenvalues are nonnegative and 
Cfc = (e]., . . . , ef) denote the eigenvectors. Therefore, the second term of the 
integrand of V t is of the form -\Y,k<d^kY,i,j<d^ifi{ u ->S u -)e k e k and is 
nonpositive due to the convexity of Q{t, •) from Assumption PO(<?). Also, 
by convexity of Q{t,-), it follows that T(x) from (12) is nonnegative and 
convex. Assumption (10) implies that the last term of the integrand of Vt is 
nonpositive and, therefore, —Vt £ A^. This implies that Q(t,St) is a local 
(.Af )-supermartingale under P. As Q{t, St) is bounded below and EQ{t, St) < 
oo, we obtain, as in the proof of the first part, that Q{t,St) is an (At)- 
supermartingale under P. □ 

Remark 2.3. 1. As seen in the proof, it suffices to assume that 



instead of (8) and (10) and, further, that So < Sq. 

2. Similar comparison results are obtained for g £ T (or g E T), where 
T £ {^st^dcxi^cx^ism^sm}. The crucial equation in the proof is (11). 
Table 1 lists sufficient conditions similar to (7)-(10) that imply Eg(ST) < 
E*g(S^) under the respective integrability conditions. 

3. For g 6 T, T £ {7" st , Fidcx , Hsm , ^icx , ^dcx , ^sm , ^cx} , similar arguments 
show that E*g(S?p) < Eg(Sr), if the corresponding inequalities of the differ- 
ential characteristics of S and S* are reversed. 
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The propagation of order property is a crucial assumption in the previ- 
ous comparison theorem. To apply the theorem, this condition needs to be 
checked in the models considered. Bergenthum and Riischendorf [3] estab- 
lished propagation of (directional) convexity for stochastic exponentials of 
Levy processes and the propagation of convexity property for some classes 
of diffusions and diffusions with jumps with zero drift coefficient and finite 
jump intensity. We establish propagation of increasing convexity for some 
classes of diffusions with jumps with nonzero drift and finite jump intensity 
and in the second part of this section, we establish the PO(.F) property for 
spatially homogeneous processes for all order-generating function classes J- 
considered in this paper. 

Let S* be a unique strong solution of the SDE 

dS; = b*(t, 5*) dt + a*(t, S*) dW t * 

(13) +<P*(t,Sl,y)(N*(dt,dy)-\*(dy)dt), 
5*0 = s 0) 

where W t * is a d-dimensional Brownian motion and N* is a Poisson random 
measure on [0,T] x R with deterministic intensity X*(dy)dt and A*(R) < oo. 
We consider the following two cases for the parameter of the jump part (j>* 
of the SDE (13). In the first case, we assume that eft* factorizes as 

(14) <l>*{t,a,y):=ip*(t,a)il)*{t,y), 

where tp* : [0, T] x M. d R + , ip*(t, •) is increasing and convex for all t £ [0, T] 
and ip* : [0, T] x R — > R rf . In the second case, we assume that d = 1 and <f>* is 
of the form 

(15) (f>*(t,s,y):= J2<pUt,s)1>i(t,v), 

i<m 

where if* : [0,T] xl-*l + , ip*(t, •) £ F icx for all t £ [0, T] and i/j* :[0,T]xl^ 
R_l_, ip*(t, •) is increasing for all t £ [0, T], for all i < m. 

Table 1 



Ordering 


Drift 


Diffusion 




Jump 




b<b* 


c — c* 




K = K* 


■^~idcx 


b<b* 


c < c* 




A < dcx A* 


■^~icx 


b<b* 


C <psd C* 




A < cx A* 


■^""ism 


b<b* 


c< c*,c" =c*" 


i < d 


A < sm A* 


F dcx 


b = b* 


c < c* 




A' < dcx A* 


F ~x 


b = b* 


C <psd C* 




A < cx A* 


F jm 


b = b* 


C < C , C = C 


,i < d 


A < sm A* 
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Remark 2.4. Under the assumption in (14), the SDE in (13) can be 
seen as a diffusion driven by a process with independent increments (PII). 
For time-independent ip* , it is a Levy-driven diffusion. In the case d = l, 
the class of functions in (15) allows the approximation of general linking 
functions (f>* and, thus, in this case, the assumption allows the treatment of 
more general diffusions with jumps. 

For to £ [0, T] and K £ N, we discretize [to,T] into K + 1 equidistant 
points ti := i T J^° + to,i £ {0, . . . , K} and denote the Euler scheme of S* 
by 

S*K,t i+1 = S*k,u + b*(U,S* K>ti )AU + a*(t h S* KA )(W t * +1 - Wft 
(16) +^(^,^,y*)iV*- J B y >*(t i ,^,y*)A*(R)At i , 

where N* is binomial with P N * = (1 - A*(K)Ati)£{ } + A*(M)A^e {1} , Y* 

has distribution ^ on (R,B), E Y ' f {-,-,¥*) := ^ //(-, ;y)\*(dy) and 

AU = t i+ i -U = Observe that W t * +1 - W t * ~ N(0, AI), where N(n, S) 

denotes the d-dimensional normal distribution with expectation \i and co- 
variance matrix £ and / £ M + (<i, M) is the identity. 

We establish propagation of increasing convexity for the propagation op- 
erator QK(t,s) of the Euler scheme that is given by 

g K (t,s):=E*(g(S* K;T )\S* Kit = s). 

We then use an approximation argument that implies propagation of in- 
creasing convexity for the propagation operator Q(t, s) also. We call this the 
approximation property. 

Assumption AP(g) (Approximation property). Let g : M. d — > R. The Eu- 
ler scheme S*^ of S* satisfies the approximation property AP(g), if for 
K — > oo, it holds true that 

G K (t,s) :=E*{g{S* K ^)\S* K)t = s)^Q{t )S ) for all t € [0,T],s£ R d . 

Similarly, the Euler scheme S 1 ^ of S* satisfies the approximation property 
AP(J r o) fo r some T§ C T if AP(g) holds true for all g £ J^q. 

The approximation property is satisfied under smoothness and linear 
growth conditions on the coefficients and smoothness and growth condi- 
tions on the test function g (cf., e.g., [14]). For suitable convergence results 
for Levy-driven diffusions, we refer to [12] and [19]. 
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Sufficient for propagation of increasing convexity for the propagation op- 
erator Gl((t, s) is <; cx -monotonicity of the corresponding transition operator 
T* (i.e., S* < icx S% implies T*5J < icx T*S%) that is of the form 

T*S* = S* + Ab*(S*) + a*(S*)W* 

(17) 

+ <j)*(S*,Y*)N* -E Y <j)*(S*,Y*)E*N*. 

Lemma 2.5 (Increasing convex monotonicity of Markov operators). Let 
S* , W* , N* , Y * be independent random variables, where S* , W* are R rf - 
valued, W* ~ N(0,AI), Aei + and N* ,Y* have values in R. Assume 
that b* : R d — > R rf and a* : R d — > M + (d, R) are increasing and convex, where 
M+(d,R) is supplied with the positive semidefinite ordering < ps d- Then the 
transition operator T* given by (17) is <[ cx -monotone if: 

1. (f)*(s,y) = (p*(s)i/j*(y), where tp* : M. d — ► R + is increasing and convex and 
'tp*:R^R d , or 

2. d=l, N*>0 and (j)*(s,y) = J2i< m V*{s)^*{y), where ipf :R^ R + are 
increasing and convex and ij)* : R — > are increasing for all i <n. 

Proof. Assume that SI and S 2 are d-dimensional random vectors that 
are independent of W*,N* and Y* and satisfy S* <; cx S%. Without loss of 
generality, by Strassen's theorem, we choose S*,S 2 such that E(S 2 \Sl) > S*. 
For / G J^icxi h follows from Jensen's inequality and independence of the 
random vectors that 

E*f{T*S 2 ) = E*E*(f{T*S 2 )\S*,W*,N*,Y*) 

> E*f(E*(S%\SZ) + AE*(b*(S%)\S*) 

(18) 

+ E* {a* (E%) | SI ) W* + E* {(f)* (S* , Y * ) \ S{ , Y * ) N* 

- E*(E Y *f(S 2 ,Y*)\Sl)E*N*). 

If (f>* factorizes into (j)*(s,y) = ip* (s)i/j* (y) as in part 1, then using the inde- 
pendence properties, the right-hand side of (18) equals 

J E*f(E*(S;\si) + AE*{b*(S^)\ Sl ) + E*(a*(S%)\si)W* 

+ E*{ l p*(S%)\s 1 ){ilj*(Y*)N* -E*(^*(Y*)N*)))P s *(d Sl ). 

Due to monotonicity and convexity of b* and a*, Jensen's inequality implies 
that E*(b*(S* 2 )\ Sl ) > b*(E*(S* 2 \ Sl )) > b*(si) and E*( Sl ) := E* (a* (S* 2 )\ Sl ) > psd 
a* (E* (S 2 )\si)) > psd a*(si). As <t*(si) and are positive semidefinite, 

it holds true that (<7*(si)) T cr*(si) < psd (S*(si)) T S*(s 1 ). A classical convex 
ordering result for normally distributed random vectors (cf. [17], Theorem 
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3.4.6) implies that a*(s\)W* < cx £* {s\)W* . In the following, we establish 
the ordering of the jump parts 

B ■= p*( Sl )(ip*(Y*)N* - E*(tp*(Y*)N*)) 

< cx E*(ip*(S2)\sx)(ip*(Y*)N* - E*(ip*(Y*)N*)) =: C. 

As ip* G ^icx, Jensen's inequality implies that $*(s\) := E* (p* (S 2 )\si) — 
<p*( Sl ) > 0. For j < d, we define Rj = ti*(si)(rf(Y*)N* - E*(^(Y*)N*)). 
Then for g 6 J 7 ^ nC 2 , it holds true that 

E*g(C) > E*g(B) + E*(Vg(B), R), 

where R = . . . , i?^), V is the gradient and (■, •) stands for the scalar prod- 
uct in R d . From R = ^rff^B it follows that E*(Vg{B),R) = P^E* (Vg{B),B). 
A characterization result of optimal couplings in [20], Theorem 1, implies 
that (B,Vg(B)) is an optimal ^-coupling, hence E*(Vg{B),B) > (E*Vg{B), 
E*B) = 0, as E*B = 0. As the convex order is generated by functions 
9 G -^icx n C 2 , it follows that B < cx C, hence B < 1CX C. 

As the increasing convex order is stable under convolutions, it follows 
from the independence of W* , Y* and iV* that 

si + Ab*(s t ) + a*{ Sl )W* + <p*{ Sl )(ip*{Y*)N* - E* (ip* (Y*)N*)) 
< icx E*(S* 2 \ Sl ) + AE*(b(S* 2 )\ Sl ) + E*(a(S* 2 )\ Sl )W* 
+ E*((p*(S 2 )\si)(ip*(Y*)N* - E*{i;*(Y*)N*)) 

and as E*f{T*S{) = jE*f( Sl + Ab*( Sl ) + a*( Sl )W* + tp*(si)(i(>*(Y*)N* - 
£*(^*(r*)jV*))P 5 i(dsi), the result follows. 

If d = 1 and <^>* has a representation of the form <f>* (s,y) = J2i< m ft ( s ) (v) 
as in part 2, then similarly to the previous case, it suffices to establish 

B := V *(sim(Y*)N* - E*0j*(Y*)N*)) 

i<m 

< cx ]T E*(<p*(S* 2 )\ Sl )(^(Y*)N* - E*W(Y*)N*)) =: C. 

i<m 

From Jensen's inequality and monotonicity of ip*, it follows that $*(s±) : = 

£*(^(S|)| Sl )-^( Sl )>0.^^^ 

and g £ J- Cx f] C 2 , Taylor's formula implies 

E*g(C)>E*g(B) + E*g'(B)R 

and it remains to verify nonnegativity of E*g'(B)R. We make use of some 
results on association of random vectors (cf. [17], Theorems 3.10.5, 3.10.7). 
As Y* and N* are independent random variables, (Y*,N*) is associated. 
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From monotonicity of tfj* > 0, it follows that **(y, n) := ip*(y)n, n > 0, is 
nondecreasing in (y,n), for all i < m and 

Z* := (Z$,..., Z*J = (*l(Y*,N*), ** m (Y*,N*)) 

is associated, as N* > 0. Thus, Z* = Z* — E*Z* is associated (cf. [17], 
Theorem 3.10.7). From nonnegativity of (p*(s%) and it follows that 

(B,R) = {J2i< m 'Pi(si)Z*,J2i<m^H s i) Z *i) is nondecreasing in Z* and, 
therefore, is associated. Thus, it holds true that E*Fi(B,R)F2(B,R) > 
E*F 1 (B,R)E*F 2 {B,R) for all nondecreasing F(.:MxR-t|, k = 1,2. As 
gGf cx nC 2 , 5' is nondecreasing and with Fi(B,R) := g'(B), F 2 (B,R) := R, 
it follows that E*g'(B)R > E*g'(B)E*R = 0. □ 

From Lemma 2.5 and Assumption AP(g), propagation of increasing con- 
vexity for some classes of diffusions with jumps follows. 

Theorem 2.6 (Propagation of increasing convexity for diffusions with 
jumps). Let g £ ^i cx and S* be a d- dimensional diffusion with jumps. Let 
Assumption AP(g) be satisfied. Ifb*(t,-) and a*(t,-) are increasing and con- 
vex for all t G [0,T], as in Lemma 2.5, and (j>* satisfies (14) or (15) (d = 1 in 
the latter case), then the propagation of increasing convexity property holds, 
that is, 

a(t,0e^icx Vie[0,T]. 

Proof. We consider the Euler approximation scheme S* K defined in 
(16) with interpolation points ti and define the corresponding transition 

operator by %.S* = S* + &*(t i ,S*)At i + a*(U,S*)W* + <f>*(ti, S*,Y)N* - 

E Y <P*{U,S*,Y)EN*, where W* = W t * +1 - W t * and P*{N* = 1) = A(R)At» 
= 1 — P*(A r * = 0). Then for to G [0, T] the Markov property implies 

g K (to,y) = E*(g(S* KtT )\S* Ktt0 = y) = E*g(T tK _ 1 • • -T^y). 

For yi,y2 G ~^- d and a G (0, 1), let Y be a Bernoulli random vector with 
distribution P Y = ot£s yi \ + (1 — aO e {j/ 2 }- Then 

oyi + (1 - a)y 2 = EY < cx F. 

Using the <j cx -monotonicity of the operator Ft for all t £ [0, T], from 
Lemma 2.5, we obtain 

G K (t ,a yi + (1 - a)y 2 ) = Gk^EY) = E*g(% K1 ■ ■ -T to EY) 

<E*g(T tK _ 1 ---T t0 Y)=g K (t ,Y). 
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Taking expectations on both sides implies convexity of Gxiior)- From 
<i cx -monotonicity of the operator If, it follows for x < y that 



GK(t ,x) = E*g(% K _ 1 • • -T to x) 

< E*g(T tK _ 1 - ■ -T to y) = g K (t ,y). 



Thus, Gic^i ") G -^"icx an d Assumption AP(g) implies that G{to, •) G ^ 



□ 



Theorems 2.2 and 2.6 imply increasing convex comparison of diffusions 
with jumps with (13) and semimartingales S ~ (b, c, K). S* has characteris- 
tics S* ~ (6*,c*,K*) of the form (3). 

Theorem 2.7 (Increasing convex comparison of semimartingales to Marko- 
vian diffusions with jumps). Let S* ~ (b*(t, S?_), a*(t,S?_)(a*(t,S^)) T , 
y*<t> (t,S t _,-)^ ^ g a d- dimensional diffusion with jumps, let S ~ (b,c,K) be a 
d- dimensional semimartingale and assume that Sq = Sq. Let g G Ficx an d 
assume that Assumption AP(g) is satisfied for S* . Let b*(t, •) and o~*(t, ■) be 
increasing and convex for all t G [0, T], as in Lemma 2.5, and let (ft* satisfy 
(14) or (15) (d= 1 in the latter case). If 



I x P-a.e., for all/el+xl^xl^R with f(t, s, •) G F cx suc/t too* tfie 
integrals exist, then 



In the sequel, we apply the comparison results of Theorem 2.2 to Levy 
processes. To that aim, we establish the propagation of order property for 
all order-generating function classes T that are considered in this paper. 
Also, in the following Section 3, we will establish comparison results for 
the finite-dimensional distributions of Levy processes, derived by a different 
method. 



We denote the transition probability by Pt t (x, B) := P*(S* — S* G B — x) 
aadi7 >t (B):=P; >t (0,B). 



Lemma 2.8 [PO(g) for spatially homogeneous processes]. Let T be a 
function class as in (1) and let g G T . If the Markovian process S* has a 



bl(cu)<b*\t,S t -(uj)) 



CtH < P sd (a*(t,S t .(u))fa*(t,S t ^Lo)) 




Eg(S T ) < E*g(S* T ). 
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Proof. Spatial homogeneity of P^ T implies that G{t,s) = J g(x + s) 

P t * T (dx). Since for / £ T, it holds that f x (-) := f(x + •) £ T for all x £ R d , 
the result follows from the stability under mixtures property of <jr. □ 

Levy processes have spatially homogeneous transition functions. There- 
fore, we obtain the following ordering result for Levy processes by combining 
Theorem 2.2 and Lemma 2.8. In this case, the differential characteristics are 
deterministic and linear in time [cf. (2)]. 

Corollary 2.9 [(Directionally) convex comparison of Levy processes]. 
Let S ~ (6, E, F) id , S* ~ (&*, £*, F*) id and S = S$ be Levy processes. Let W 
and W* be defined as in (5). 

1. Let g £ ^"idcx be bounded below or g £ ^idcx and f or v := F x \, v* := 
F* x \, assume that (\W\ * v)t, (\W*\ * u *)t £ ^\ oc - Then the conditions 

b<b*, Z ij <T,* ij Vi,j<d and F < dcx F* 

imply that Eg(ST) < E*g(S^). 

2. If g £ J-i cx , then the conditions 

b<b*, E< psd £* and F < cx F* 
imply that Eg(Sx) < E*g(S^). 

Remark 2.10. According to Remark 2.3, similar statements hold true 
for the orders generated by F st , ^cx^dcx^icx^idcxj-^sm and :F ism . 

In particular, Corollary 2.9 and Remark 2.10 imply the following, inter- 
esting ordering results for multivariate normal random variables which are 
established in a different way in [17], Chapter 3.13: 

Corollary 2.11 (Ordering of normal random vectors). Let N ~ iV(/i,E), 
N* ~ (/**,£*). 

1. 7/ /i < /i* and £ = £*, i/ien A < st iV*. 

2. //> = ^* and £ < psd E*, i/ien A < cx A*. 

3. 7/ /U = n* and E^' < £**•? Vt, j < d, then N < dcx A*. 

4. I//i< ii* and E < psd S*, then N < icx TV*. 

5. If n < n* and E^ < £**•? Vt,j < d, then A < idcx A*. 

3. Comparison results for Levy processes. In this section, we investigate 
a different approach that yields comparison results for finite-dimensional 
distributions of Levy processes. For this approach, we first establish an or- 
dering result for time-homogeneous Markov processes based on a mono- 
tone separating transition kernel. We consider the finite-dimensional or- 
dering (Si) < T (Sf ) based on the function classes := {/ := (R d ) m 
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M : f(s\, . . . , Si-i, •, Si+i, . . . , s m ) G F, Si G M d , i < m}, m, d G N, that are com- 
ponentwise in J 7 . A d-dimensional process is said to have smaller finite- 
dimensional distributions with respect to the product ordering induced by 
T than a d-dimensional process S (2) , (S\) <t (S 4 2 ), if for every m G N and all 
< t x < ■ ■ ■ < t m < T, it holds true that Eg(S$ ,...,SjP)< Eg(Sg> , . . . , s£ } ) 
for all 9 GfW. 

For time-homogeneous Markov processes, the existence of a <^-monotone 
transition kernel that separates the transition kernels of is sufficient to 
establish ordering of the finite-dimensional distributions. This is stated in 
the following separation result which holds true for general integral orders 
induced by a function class T\ 

Proposition 3.1 (Ordering of finite-dimensional distributions of Mar- 
kov processes) . Two time- homogeneous Markov processes (<% )te[t 1; T] an< ^ 
)te[ti,T] w ith transition kernels Q^p and Q| 2) satisfy 

(S 4 (1) )<^ (2) ) 

if S^ 1 < jr and if a family (Qt) of -monotone transition kernels exists 
such that 

qP(x, •) <jr Qt{x, •) <jr QP (x, •) for all x and all t G [ii, T]. 

Proof. The proof uses arguments similar to those used in the proof 
of Theorem 5.2.15 in [17]. We consider the case of dimension m = 2. For 
/ G J 7 ^, <jF-monotonicity of Qt 2 -t 1 {si, •) is equivalent to g(s\) := 
f f(si, S2)Qt 2 -t 1 (si, ds2) G T. As /(si,-) G the ordering of the transi- 
tion kernels implies for f^(si) := / f(s\, S2)Q[^ tl ( s i; ^ s 2) that f^(si) < 
g{si) < / (2) (si). Hence, s£ ] < T S$ implies that 

Ef(S^\s^)=J /«( Sl )P<Vi)< / g(si)P s ^ (d Sl ) 

< I g( Sl )P S 'i (dsi) < | /( 2 )( Sl )P< ) (^ 1 )= J B/(4 2) ,5; 2) ). 
The result for m> 2 follows by induction. □ 

In order to apply Proposition 3.1 to Levy processes, we consider, in the 
first step, the particular case of compound Poisson processes. We make use 
of a representation of compound Poisson processes as random sum processes 
with a Poisson number of summands, which allows comparison results with 
respect to the order-generating functions T from (1) by a natural coupling 
argument. In the second step, we extend the comparison result to pure jump 
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Levy processes with infinite Levy measure. We truncate the corresponding 
Levy measures around the origin and give ordering conditions that imply 
ordering with respect to the order-generating functions in (1) of the induced 
compound Poisson processes. By an approximation argument, this implies 
ordering for the limit process. 

For finite measures i = 1, 2, we write < t to denote that 

f f{x)M ( - 1 \dx) < J f(x)M<V(dx) for all / G JF. 

Lemma 3.2 (Ordering of compound Poisson processes). Let ~ {W' (0), 
0,F®)q } let E\sf' \ < oo, and assume that the Levy measures have the 
same finite total mass A. If F^ <p F^ holds true for: 

1. ^G{^cx,^dcx,^"sm} andbW(0)=b^(0) 
or for 

2. T G {F s t,J-icx,3~idcxi ^ism} and 6«(0) < b^(0), 
then (5 t (1) )<^(5 t (2) ). 

Proof. Let*e[0,T]. 

1. Let T G {^cxj^dcxi-^sm}- Then are representable as random sum 
processes 

N t 

S f) = 6«(0)t + £Af , * = 1,2, 
i=i 

where (Xj) are i.i.d. with distribution R^\dx) = jF^(dx) and are inde- 
pendent of the Poisson process N with intensity A. Thus, S^ 1 ' and are 
naturally coupled by the same Poisson process N. As \\F®\\ =A, condition 
< r F^ implies that <jr x] 2) Vj. From the stability of < T with 

respect to convolutions and mixtures, the coupling of and implies, 
for / G T, that, 

^/(E^r) =^ t ^/(E^ 1) ) <E N *Ef(±X^ 

= E^Ef(±xf^=Ef(j:xf^ 

where E Nt denotes the expectation with respect to the distribution of Nt. 
Also, from convolution stability, it follows that sj: ^ <jr . 
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2. Let T G {^"st, -^icx, -^idcx, -^"ism}- Similarly to the first part of the proof, 
it follows that Y^=i Xj 1 ^ <t 2~Z j^i ■ As an functions in T are increasing 
and (0) - &( 1} (0) > 0, it follows that 

Nt N t 

S« = feW (Q)t + *f ] <T b {1) (0)t + J2 Xf ] 

N, 



^{2) _ c (2) 

3. It remains to prove that the finite-dimensional distributions are also 

ordered. By Proposition 3.1, we must establish that q\ (a;, ■) = P{S^ G 

• |So = x) is <jp-monotone. To that end, it suffices to prove that for / G J-, 
the function 

f Q m(x) := J f(y)Qi 1 \x,dy) 
belongs to T . From spatial homogeneity, it follows that f n (i) is of the form 



and a G (0, 1), this implies that 



/ Q (i)(z) =Jf(x + v)Qr(dy), where Q t w (dy) :=Q^(0,dy). For T = F a 



f nW (ax+(l-a)z)= I f(a(x + y) + (l-a)(z + y))Qi 1) (dy) 



< «/ Q (i) (x) + (1 - a)/ Q (i) (^), 

thus / (i) G JP"cx- Similarly, it is established for T G {J^st^icx^idcx^ism, 
•Fdcxi-Fsm} that / (i) G T holds true. □ 

We now turn to ordering results on Levy processes with infinite Levy 
measure. We truncate the corresponding Levy measures around zero, use 
ordering results from Lemma 3.2 and then obtain orderings for the limit 
processes by weak convergence. 

Let FW be Levy measures with infinite total mass. For sequences T 
and £n I 0, we define the truncated Levy measures F« by 

(19) F®(dx):=l w _« (x)F®(dx). 

For (0,0,4° )o, we introduce the modified truncated Levy measures 

F® defined as 

(20) F«(dz) := F®(dx) + (\\F^\\ - ||F« || )H_ C|| ^3-^ „ ^ „ } 5 {0 > (^) - 
Then ~ (0,0,F"J°) and pf^ = HF^I]. 
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Theorem 3.3 (Convex- type comparison of Levy processes). Let T G 
{^cxj^dcxj-^sm}- Let sf\ i = 1,2, be Levy processes with E\s[ \ < oo and 
~ (ESi\o, F^)id, where the Levy measures have infinite total 
mass. Let e„ f and e$ [ be sequences such that for Fn given in (19), 
it holds true that J xFn\dx) = J xFn\dx). 

If ES^ = ES^ and if the modified truncated Levy measures F^ defined 
in (20) are ordered by 

(21) F^<?F ( ? VnGN, 

then (SP)<F (SP). 

Proof. Let £$ ~ (ESP,0,Ff) id = (& W (0)), 0,7f)o, where 6^(0) := 
ESf ] - JxF^(dx). Then 6^(0) = &i 2) (0) and due to Lemma 3.2, the or- 
dering on the modified truncated Levy measures (21) implies that < T 

(S$). As ES® t = ESf \ it remains to show that (S® t ) " (S^) (see [17], 
Lemma 3.4.5, Theorems 3.4.6, 3.12.7 and 3.3.12). A sufficient condition for 
functional weak convergence is convergence of the differential characteris- 
tics, as specified in [13], Corollary VII. 3. 6. For the drift component, we have 
bn\h) =b^\h), for all n. For the modified Gaussian characteristics, it fol- 
lows that 

c« M = J h k (x)h\x)F^(dx) 

= J h k (x)h\x)t (£ ^_^ )c {x)F^{dx) 

^J^ M F^ dx) ^\ 

due to the Lebesgue theorem. The appropriate convergence of the Levy mea- 
sures is with respect to g G C2(M d ) := {/ : M d — > M, / is bounded, continuous 
and has value around 0} . It is obvious that there exists an JVgN such 
that 

F®{g)= J g(x)t { _ £ ^^ ))c (x)F^(dx) 
= J g(x)F®{dx) 

for all n > N . Therefore, Jacod and Shiryaev ([13], Corollary VII. 3. 6) implies 
functional convergence 

S-W JL> thus weak convergence of the finite- 
dimensional distributions, and the result follows. □ 
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From Theorem 3.3, it follows that if a Levy process has higher activity 
in every jump height than some further Levy process then < 

cx 

(S t (2) ). 

Corollary 3.4 (Domination criterion). Let 5 t W , S® ~ (-ES^, 0, F«) id) 

i = 1,2, 6e one- dimensional Levy processes with E\S{ \ < oo and assume 
that Ft 1 ' have densities /W s^c/i i/icrf |x|FW((ix) = oo /or A = (—1,0) 

and A = (0, 1). // ES[ 1] = ES{ 2) and 

(22) 0</ {1) (x) </ {2) (x) VxGM, 
^er» (5 t (1) )< cx (,S t (2) ). 

Proof. Let g$ T 0, i = 1,2, be sequences with e ( n 1} < e ( „ 2) , for all n £ N. 
Due to the assumption J A \x\F^\dx) = oo for A = (—1,0) and A = (0, 1), 

there are sequences En I 0, i = 1,2, such that := / xF n l \dx) = 0. For 
ij\ 

the finite measures F n given in (20), define the Levy distribution functions 

F n (x) := F n ' ((— oo,x]). We establish that for all n £ N, there is an x n G R 
with 

(23) ^(x) <F {2 \x) Vx<x n and (x) > F {2) (x) Vx>x n . 

As En = 0, the cut criterion by Karlin and Novikov (cf. [17], Theorem 

1.5.17) implies F n < cx F^ 2 for all n £ N and the result follows from Theo- 
rem 3.3. 

From f^(x) < / (2) (x) and <e ( n 2) , it follows that F%\x) <F {2 \x) 
for all x < and it remains to prove that the cutting point x n is in the 
positive half-axis including the origin. The Levy distribution functions F n 
depend on e n and \\F n \\. The following three cases are possible: 

1. E n S £n , \\r n || ^ ||i*n ||, 
9 (1) (2) ,|p(2)|, ,|p(l)|| 
I. E n S £n , || -Tn || S ||-Fn ||, 

3 E {1) >E {2) 

Observe that e n 1} > e (2) implies ||F,1 2) || > \\F^ ||. 

We establish (23) for the first case; the other cases are similar. From 

ll-^n^ll > ll-^n^ll> it follows that the modified truncated Levy measures F^ 
are of the form F^ = F n 2 ^ and F^\dx) = fP (dx) + (||F,1 2) || - \\F r ( P \\)5 {0} (dx) 
and that one of the following two cases holds true: 

F^iO) <F®(0) or Fi 1) (0)>Fi 2) (0). 
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The first case implies that (x) < (x) for all x G [0, e^]. From \\F® || = 

A n and f^(x) < / (2) (x) for all x G R, it follows that F%\x) > F%\x) for 

all x >En^- As is constant on [e^e^] and F^ is strongly mono- 
tone in this interval, it follows that (23) holds true for one and only one 

%n G [En ^ , Sn ] • 

In the second case, it follows from F%\o) > F%\o) and < £ ( n ] that 

^(s) > F { *\x) for all x G [0,£^ 2) ]. The condition f^(x) < f^{x) for all 

x G R, together with \\F^ \\ = X n , implies that F^\x) > F%\x) for x > e { n ] '. 
Therefore, (23) holds true for x n = 0. □ 

Corollary 3.4 can be applied to compare one-dimensional normal inverse 
Gaussian processes in the shape and scaling parameters a and 5. The Levy 
density of an NIG=NIG(a, (3, 5, pi)-distributed random variable S is given by 

SaK^alxDeP* 



TT X 



where K\ denotes the modified Bessel function of the third kind with index 
1, a > 0, < 1/31 < a, and S > and S has expectation ES = a H — ^= . In 

a 2 —/3 2 

the following example, where we give an ordering in the parameters a and 5, 
respectively, we need equality of the expectations ES^ , which is obtained 
by choosing //W appropriately. 

Example 3.5 (Convex ordering of one-dimensional NIG processes). 

1. Ordering in a. Let S®,i= 1,2, be NIG processes with S® ~ (-E-Sf^O, 
f( l \x)dx), where for / Qj/ g i( 5 given in (24), we define f^\x) := f a (i) ps( x )- 

If a« < q( 2 ) and £S{ 1} = £Sj 2) , then (sf } ) < cx (S 4 (1) ). 

2. Ordering in <5. Let 5W,t = 1,2, be NIG processes with S® ~ (-ES^.O, 
f^(x)dx), where for / Qj/ j )<5 given in (24), we define /W(:c) := f a ,/3,5(^( x )- 

If S« < 5( 2 ) and = £Sj 2) , then (S t (1) ) < cx (S t (2) ). 

Proof. 1. For fixed x > 0, consider g(a) := f a ,(3,s( x )- Then </(a) = 

-^a-K (ax) < 0, that is, f^(x) > f^{x) for all x G R+. For fixed x < 0, 

we similarly obtain </(a) = -^-aKo(-ra) < 0, thus f^\x) > f( 2 \x), for 

all x G R_. Corollary 3.4 implies that (fif >) < cx {s\ 1] ). 

2. It is obvious that <5 (1) < <j( 2 ) implies f w {x) < / (2) (x), for all x G R and, 

therefore, by Corollary 3.4, it follows that (sjp) < cx (S t {2) ). □ 

The following result for increasing- and increasing convex- type orders is 
similar to Theorem 3.3 and Corollary 3.4: 
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Theorem 3.6 (Increasing convex-type comparison of Levy processes). 
Let T G {^stj-^icxj^idcxj^ism} and assume that sj: l \i = 1,2, are Levy pro- 
cesses with E\s[ i] \ < oo and S® ~ (£<sf ] , 0, F«) id , w/iere ||FW|| = oo. As- 
sume that ES^ < ES^ and let e$ f an d £n i be sequences such that 
for Fn given in (19), it holds true that < / xFn\dx) — J xFn\dx) < 
ES{ 2) -ES[ 1] . 

If the modified truncated Levy measures defined in (20) are ordered 

as 

then (S[ l) )<r {s[ 2) ). 

Proof. For 6^(0) := ESf -JxF®(dx) and S$ ~ (6^(0), 0,F?) , it 
follows from Lemma 3.2 that (S$) <f (5$). As ES^ = 6^ (id) = £5?, 
it remains to prove functional weak convergence. This is similar to the proof 
of Theorem 3.3. □ 

If the processes have paths of finite variation, then there is a similar 
comparison result that does not require the first moment condition < 
/ xFn^ (dx) — J xFn^ {dx) < ES^ — ES^ . In this case, we postulate an 
ordering condition on the drift terms (0) < b^ (0) of £W . The ftW (0) 
exist, due to the fact that Levy measures of processes with paths of finite 
variation integrate \x\ around the origin. 

Theorem 3.7 (Increasing convex- type comparison of Levy processes with 
paths of finite variation). Let T € {^"stj 3~icx-> -^"idcxj -^ism} o,nd assume that 
sf\i = 1,2, are Levy processes with paths of finite variation. Assume that 
E\sf\ <oo and S® ~ (6 (i) (0), 0, F«) , where \\F®\\=oo. Let T and 
£n I be sequences such that for F$p given in (19), it holds true that 
fxF^\dx)<fxF^\dx). 

If b^ (0) < b^ (0) and the modified truncated Levy measures F^ defined 
in (20) are ordered as 

F^<fF {2) VnGN, 

then (Si 1] )<^(si 2) ). 

PROOF. For S$ ~ (&W(0), 0,^)0, it follows from Lemma 3.2 that 
(Sg) <t (S$). As ES% = &£>(id) = &«(0)+ frf(dx) - &«(id) = £Sf\ 
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it remains to prove functional weak convergence. This is similar to the proof 
of Theorem 3.3. □ 

For the one-dimensional Theorem 3.7, a sufficient condition for the or- 
dering with respect to J- S t in terms of Levy densities. If the Levy density of 
one of the processes is dominated by the Levy density of the other on the 
negative half-axis and the domination is reversed on the positive half-axis, 
then the first process is smaller than the second one with respect to the 
usual stochastic ordering. 

(i) 

Corollary 3.8. Let SI ,i = 1,2, be one- dimensional Levy processes 

with finite variation, E\s[ \ < oo and S® ~ (b^ (0),0, F®)q. Assume that 
FW have densities /W that are monotonically increasing to infinity as x 
tends to zero. 

I/6 (1) (0) <& (2) (0), 

f {1 \x) > / (2) (x) VxGR_ and f (1 \x)<f {2 \x) VxGlR + , 
then (S^)< st (Sj; 2) ). 

Proof. Let e n J. and let the S6CJU6I1C6S £_ /I ) T 0, ei ] | in (19) be given 
by := £ni £ n '■= —tn, i = l,2. Then for En := J xFn\dx), it follows that 
E n < E n , due to the pointwise ordering condition on the Levy densities 
/W. Similarly to Corollary 3.4, we establish that 

(25) F^\x)>F ( Z\x) Vxel. 

This implies that F^ 1 < s t F^ , for all n G N and the result follows from 
Theorem 3.7. 

From > f (2 \x), for all x < and = -e n it follows that F^ (x) > 

(2) 

F n (x) for all x < and it only remains to prove the domination criterion 
(25) on the positive half-axis including the origin. The modified truncated 
Levy measures depend on the sign of \\F^\\ - \\fP\\. If ||Fi 2) || > IIF^H, 
then = Fi 2) and F%\dx) = F n 1] \dx) + (||^ 2) || - \\F n l) \\)8 {Q} (dx). This 
implies that F {1) (x) >F {2 \x), x G [0,e n ]. From / (1) (x) < / (2) (x), for all 
x G R + , it then follows that F^\x) > F^\x) also for x > e n and (25) fol- 
lows. If || < ll-^n^ll then a similar consideration implies (25). □ 

Due to the assumption of finite variation, the previous comparison crite- 
rion is not applicable to two NIG processes with different skewness param- 
eters 0^' < /?( 2 ) , although the domination criterion on the Levy densities of 
Corollary 3.8 is satisfied. 
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4. Ordering result for mixing-type distributions. For some cases of in- 
terest, it is possible to obtain comparison results by using mixing- type rep- 
resentations. We apply this approach to the class of generalized hyperbolic 
(GH) distributions. Further mixing-type representations of particular inter- 
est for financial mathematical models hold for multivariate ^distributions 
and elliptically contoured distributions (see [4]). GH distributions are vari- 
ance mixtures of multivariate normal distributions which have a general- 
ized inverse Gaussian distribution as mixing distribution. For £ WL d , 
A« £ M(d,R) with det(A^) = 1 for i = 1,2 and N® ~ M{0, AW), we con- 
sider the d-dimensional random variable 

(26) S {i) = fi {i) + X« A® /?« + iV (i) , 

where X" are generalized inverse Gaussian random variables with densities 

(27) c* GIG(M7)( *) := (I)"^^^-^)^/^)!^^), 

where 5 > 0, a 2 > (3A(3 T and 7 = \J a 2 — [3A[3 T . Then 5 is generalized hy- 
perbolic distributed with parameters d, A, a, j3,S, \jl and covariance matrix A 
and we write GH(d, A, a, (3, 6, fi) (cf. [1]). 

The following lemma states a comparison result for GIG distributions 
with respect to the likelihood ratio order <i r , if the parameters A, 5 and 7 
are ordered: 

Lemma 4.1 (Likelihood ratio ordering of GIG random variables). Let 
X^> be GIG distributed with density c^gig(aW 5« y(.i)){ x )- 
If A« < A< 2 > , <5« < 5( 2 ) and 7 « > 7 ( 2 ) , then X« < lr . 

Proof. We consider the likelihood ratio 

, s d GIG(A« ,5(i) n W)( x ) 

9{x) :=-j r -r 

«GIG(A( 2 ),5( 2 ),t<( 2 ))W 

= ^ (1) - A(2> e( 1 /2)(((^ 2 -^ 1 ) 2 )/,) + (7' 2 ) 2 -7( 1 ) 2 )x) lM+(2;)) 

with K > 0. The first derivative of g is of the form g'{x) = K\(x)(\^ — 

\ (2) )+K 2 (x)( 8W2 -f )2 + ( 7 ( 2 ) 2 - 7 W 2 )), where K l (x),K 2 (x) >0. Thus, it 
follows from the orderings on the parameters that g'(x) < 0. □ 

In the following theorem we give an increasing convex comparison result 
of multivariate GH distributions with mixing- type representation (26). We 
consider the following three cases for 0(0 and A»: 

(28) 0</3 (1) </3 (2) , A ( = /, 
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(29) 00=0, A«< psd A( 2 ), 

0<{P)<fP\ AW< psd A( 2 ), 

(30) 

< A^ < A^ Vt,j<d. 

Theorem 4.2 (Increasing convex comparison of GH distributions). Let 

S« be GH(d,AW,aW, 0(0, jW^W) distributed. If 

A (D< A (2) ; < ( j(2) ) a (l)> a (2) s ^(D<^) 

and one o/ i/ie cases (28)-(30) ZioWs true /or 0(0 and AW , i/ien S^ 1 ) <i cx 

Proof. First we prove that the conditions on the parameters of 
imply that X^ < st X( 2 ). For 7W = ^(O 2 - 0(OA(O(0(O)7\ i t follows from 
the conditions on aW, AW and 0(0 i n a n three cases (28)-(30) that 7W > 
7( 2 ). Due to Lemma 4.1, the ordering conditions on A(0 and JW imply that 
X« < b X( 2 ) and this yields X® < st X( 2 ). 

1. Let 0(0 and A« satisfy (28). For / £ ;F icx , we define := £/0uW + 

pWx + y/xN), with iV = iV(0, f = 1,2. Let < < x (2) . As p^x + y/xN ~ 
M((3^x,xl), 0( 1 M 1 ) <fi0) x W and sWj < psd x (2) I, it follows from Corol- 
lary 2.11 that g is increasing. Then XW < st A"( 2 ) implies that 

Ef(fi^ + 0«A:( 1) + \ZxWN) = Eg(xW) < Eg(X®) 

= Ef(jF> + 0«X( 2 ) + \fW)N). 

As / is increasing and X ( 2 ) is nonnegative, it follows from /Z 1 ) < /i( 2 ) and 
0(1) < 0(2) t hat 

jE/ ( /U (i) + 0(i)x( 2 ) + Vx(2Tiv) 

= E x(2) Ef(^ + ^(i) x (2) + y/^)N) 

< E xm Ef(^ + 0( 2 )x( 2 ) + Jx^N) = Ef(SW). 

2. Let (i) and A(0 satisfy (29) and assume that / € ^"i cx and > 0. 
From AW < psd A( 2 \ it follows that x^A^ < psd x^A^ and, therefore, 

Corollary 2.11.4 implies Ef(S^) < £/(u {2) + vfW N (2) ). For z > 0, let 
g(z) := i?/(/i( 2 ) + zA r ( 2 )). As A( 2 ) is positive semidefinite, it follows for 
< z (2) t h a t 2 (!)A( 2 ) < psd z( 2 )A( 2 ) and Corollary 2.11.4 implies that g 
is increasing. From X^ < st A"( 2 ), it follows that 

Ef{^ + \fx^N^) = Eg(^5(W) < Eg{\fx&) 
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3. Let and A® satisfy (30). The condition A$ > 0, £, j < d, implies 



that A® pW T < fiP) A® P& 1 for < /3« < /?( 2 ). The result follows sim- 
ilarly to the previous ones. □ 

In the case A = — i, is normally inverse Gaussian distributed. NIG 
distributed random variables are stable under convolutions: 

NIG(d, a, P, 5, n, A; t) = NIG(d, a, P, tS, tfi, A; 1). 

Therefore, Theorem 4.2 also implies increasing convex comparison of the of 
NIG processes with mixing-type representation 

(31) := n®t + xPa^P® + y/x?>N®, 
where xf ] ~ GIG(-|,ti5W,7W) for any time i > 0. 

Corollary 4.3 (Increasing convex comparison of NIG processes). Let 
Sf ] be NIG(d, aW , 00 , $0 , ^0 , A « ; t) processes. If 

and one o/ £/te cases (28)-(30) /ic-Zds irne for /?0 and AW, i/ien (5 t ) <i cx 
(5 t (2) ). 

Proof. As in Theorem 4.2, it follows from the ordering conditions that 
St <icx for all t € [0, T]. The result on finite-dimensional distributions 
follows from Proposition 3.1. □ 

APPENDIX 

Proof of Lemma 2.1. Let t e (0,T] and assume that 

(32) / |Ag(n,^_(a;*),n)|^(^;dn,dn)G< c . 

Ito's lemma implies that Q{t,S*[) is a semimartingale with evolution 
g(t,$)=g(p,SZ) + M t + I (Ag)(u,S:_,y) f i*(du,dy) 

J[0,t]xR d 

+ / |D^(n,,S:_) + ^D^(n, > S:.)6"(n,S:_) 

+ D 2 ^(n,5;„)c^(n, < S:„))dn, 
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where M t :=J2i<dI\ot]^- > i^( u ^u~)dM*' 1 is a one-dimensional local (A%)- 
martingale under P* and M* % denotes the ith component of the martingale 
part of S* . As (32) holds, there is a local martingale Mt such that 

Git, S* t ) =G(0,S* ) + M t + M t 

+ / \D t G(u,SU+J2 B £( u > S u-)b n (u,S:_) 

+^ D ^(«,tK ij (^:-) 

i,3<d 

+ f (AG)(u,S:_,y)K* u (S:_,dy)\du 

=:G(0,S* ) + M t + M t + V t , 

where we denote the Lebesgue integral by Vt. As G(t, Sf), M t and M t are 
local (^)-martingales under P*, Vt is a local martingale starting at zero, 
and as Vt is of finite variation, it follows that V = 0. Therefore, we obtain 

D t G(t, SI) + D ^(*> #_)&**(*, St-) 



+ i J2 VlG(t,S:_)c* l >(t,Sn+ I (AG)(t,S:_,y)KZ(Sl,dy) = 
and (4) follows. 

It remains to prove that Q(t,-) £ J- Cx implies the integrability condition 
(32). As Q(t,Sf) is a local martingale under Q* , it is a special semimartin- 
gale. The process 

/' D t g(u,S* u _)du + ( YPiG(u,&LJ)dA« 

is predictable and of finite variation and is therefore in -4i oc . Due to a repre- 
sentation result for special semimartingales (see [13], Proposition 1.4.23), it 
follows that (W* *fJ,*)t G A\ oc . As convexity of G(t, •) implies AQ(t,S^_,y) > 
0, (32) follows. □ 
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